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ABSTRACT 

/ 

A matrix method using flexibility influence coefficients was^''^'"^ 
developed for obtaining the free-free bending and torsional mode 
shapes, slopes of mode shapes, and natural frequencies of space 
vehicles. The effects of rotary inertia and shear flexibility are 
included. 

The mode shapes and natural frequencies were determined for a 
typical space vehicle and compared with those obtained from a modified 
Stodola method. Twenty mass points were used for the influence coef- 
ficient analysis and 201 for the Stodola method. 

This report shows that the influence coefficient method will 
obtain satisfactory mode shapes and frequencies in comparison to a 
Stodola method. 
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SUMMARY 

Mode shapes and natural frequencies of a uniform and nonuniform 
beam obtained by a flexibility influence coefficient method were compared 
with a modified Stodola method. Twenty mass points were used with the 
influence coefficient method, while 201 were used with the Stodola method 
for the nonuniform beam analysis. For the uniform beam, 10, 15, and 20 
mass points were used, while 145 were used in the Stodola analysis. The 
results indicate that the flexibility influence coefficient method yields 
reasonably accurate mode shapes and frequencies in bending and torsion 
using only twenty mass points. Maximum variation in frequencies between 
the two methods used was only about 2 percent for the nonuniform beam and 
3 percent for the uniform beam. Mode shapes compared very closely with 
the exception of the 4th modes of the free-free bending and torsion for 
the nonuniform beam. 


I. INTRODUCTION 

During the past few years the vibration of various structures 
and their components has become increasingly important to scientific 
personnel in many fields. Practically any structure which is subjected 
to shock or repeated loads experiences vibrations. These vibrations 
result, in many cases, in structural fatigue, due to repeated stress 
reversals, or violent structural failure due to a resonant condition. 

In the space field, various problems arise in the design of con- 
trol systems because of the elasticity of the structure. Insulation 
of sensitive instruments against shock and vibration is a problem which 
must be considered. Also, acoustical problems arise due to the high 
energy level of the sound waves emitted by the rocket motors. 
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This report is concerned specifically with the vibration of non- 
uniform beams, a problem which is analogous to the structural vibration 
of a space vehicle airframe. A matrix method using flexibility influence 
coefficients to determine the mode shapes, slopes of mode shapes, and 
natural frequencies, both torsional and bending, for a uniform and non- 
uniform single-beam structure is presented. This analysis includes the 
effects of rotary inertia and shear flexibility. 

The author expresses his appreciation to Mr. C. R. Wells of 
Chrysler Corporation Space Division for the many helpful suggestions 
in the preparation of this report. 


II. DESCRIPTION 

The total linear or angular deflection of any point on a beam can 
be expressed as the sum of the deflections at that point produced by 
the individual applied forces and torques. This is the principle of 
superposition which will be used in writing the deflections and slopes 
of a vibrating beam. The general equation for the displacements or 
rotations of points on a beam can be written in the following form: 


n 

qi = ^ c. j q. (i = 1, 2, 3,. . . n) (1) 

j=l 

where q^ 1 s are the generalized coordinates, deflection and rotation, 
Cij’s are flexibility influence coefficients, and Qj’s are generalized 
forces or torques. The flexibility coefficients can be determined by 
subdividing a beam into n parts, assuming the mass of each element to 
be concentrated at the center of the element, applying a unit force and 
moment separately at each point, and then determining the deflection and 
slope at each point on the beam for each loading condition. Influence 
coefficients of this type can be thought of as the reciprocal of the 
spring constants for each mass point. The generalized forces Qj are 
the inertia forces miYiw 2 and the inertia torques IjY^co 2 . 

There are seven types of flexibility influence coefficients 
associated with bending vibration problems. Equations (2), (3), and 
(4) illustrate their relationship to the total deflection, slope, and 
bending slope of a beam. 

From equations (1) and (2), the total deflection and slope can be 
obtained for the i fc h mass point along a vibrating beam. 
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y ti - -^K^j + °n) “i y ti + + Y n 


r— i / 


+ ^^ + <Ki y m 


( 2 ) 


v 1 / 


Y ii ‘ ^L(, 9 ” + 9 «)’i Y ™ + “"I 9 ^ J « Y ii + **i -«l- 


I Y' , 

x r> • j ‘t» • 4 


(3) 


Equations (2) and (3) have 3 unknowns, Y^,^, Yi,^, and Y^; therefore, 
an equation for Yfci must be written before a solution is possible. 




FB 

0. . m. Y 
ij i Ti 





q^® i Y* 
lj Ri Bi' 


(4) 


Two different mass moments of inertia are used in equations (2), 

.(3), and (4): l£± and I R ^. Igi is the inertia due to the length of 

the section and I R ^ *- 8 the inertia due to the radius or diameter. In 
equations (2), (3), and (4) one sees that I^ always occurs with Y^, 
while I R £ occurs with Y^. The reason for this is that the length of an 
element of beam rotates under both bending and shear loads; therefore, 

Igj_ is always accompanied by Yi;i + Yj^ or Y^^. In the other case, the 
diameter or radius of an element rotates only when subjected to a bending 
load. A shearing load causes sliding of adjacent planes in the vertical 
direction, but does not produce any rotation of the diameter of the 
element with respect to the vertical; therefore, I R £ is associated only 
with Yg^. 

Since the center of gravity of each element may not coincide with 
its geometric center, it is necessary to add additional forces and 
torques due to this unbalance as follows: 


(1) A force: co 2 nu 

(2) A moment: co 2 nu 

(3) Another moment: co 2 nu S? Y,^. 

Equations (2), (3), and (4) now read: 



JFS 


+ C“ m Y . + S Y' + uf 


ij ij J “i V Ti i ‘Ti 




+ <%) ( J *i + “i S f>' 


+ -V L Ri Y Bi 


(2') 


■ “"K 6 !^ e ij) m i( Y Ti + S i Y ii 

S i Y Ti + M2 I e “ 


+ co 2 > 0*°? ( J„, + m. s? ) Y' 


I . Y' 
Ri Bx 


ij V -ei i V Ti 

(3') 


= co 


m i ( Y Ti + S i Y ii) + “"I ^ ( J *i + "i S i) Y Xi 


t- co" 




(4*) 
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Equations (2* ), 


(3') and (4') 


can be written in matrix form as 


follows: 



T 



I 


l 1 J 



Since flexibility influence coefficients are more easily obtained 
for a cantilever beam than beams with other end conditions, it will be 
assumed that the above coefficients have been determined for this case. 


Equation (5) will be written as 




( 6 ) 


and upon iteration the mode shapes, slopes of mode shapes, slopes of 
the bending mode shape, and natural bending frequencies are obtained. 

To obtain the free- free frequencies, modes shapes, etc., the clamped 
end of the beam must be allowed to translate and rotate as shown below. 




The cantilever deflections and slopes and 
written in terms of the new variables Y^, x i> ®o» 


can now be 
and Y 0 . 


Y 


Ti 



X. 0 
1 o 


( 7 : 




( 8 : 



(9: 


where Y^ is the deflection of the free-free beam and Y^, is the slope 
Matrix equation (5) now reads 



\ 
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O 1 o 
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1 
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1! 
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Y'. 

, L 

Y Bi F 

CD 

o 
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Y» 
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1 > 



V > 
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(10 


Writing equation (10) in three separate equations, 



(11a 


(lib 
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where 


ce 


Jr i 

1 


is the first partitioned row of the influence coefficient 


matrix | C0 j the second row and 
*2 


I- ° Jr s 

moment of inertia matrix is 
designated 


ce 


-1r 3 


the third row. The mass and 


ml 


M 


. The unknowns Y., Y’ and Y* . are 
i l Br 


{■}• 


Introducing the boundary conditions that the shear and bending 
moment are zero at the ends of the beam, we obtain 


X m i ( Y lF + S 1 Y i F ) - 0 (12) 

I m i ( x i + s i) (> + s i n P ) + 1 ^ n F + Z ^ Y h F - 0 < i3 > 


where 


1 cg * hi - m i s i • 


In matrix form, equations (12) and (13) read 



( 12 *) 


( 13 * )- 



Solving equations (11a), (lib), and (11c) for Y. , Y'. 
substituting in (12') and (13') yields the following: 


, and Yj 



and 



For simplification, the following substitutions are made: 






“ 


- 




- 


- 

A 

= 

m i 


ce 


ml 

+ 

m. S . 


ce 


ml 

- 




- 

R i 

- 


i i 


_ 

r 2 

- 






“ 






- 


- 

B 

— 

m. X. 

1 X 


ce 

R i 

ml 

+ 

m. S . 
l i 


ce 

R i 

ml 


D 


ce 


T ml 1 + 


m. X. S. 


[ ce 


ml 



9 


Equations (14) and (15) may be written with the above 
as follows: 


Y 0 ^ m i + ®o ( X i + s p + tj2 A u J- = 0 . 

y q £« t <x t + s.) + e 0 [^- t Of + 2x t Sl ) + + £ 



Making further substitutions: 


M o = ^“i 

x o - a»i (Xf + 2x £ s t ) 



S c = 2 * 4 . (X t + S t ) 


L o = S < J *i + J Ri>* 


Equations (14') and (15') are written in the following form: 


M o Y o + S o 6 o + w 2 


u j- =0. 


S Y + (J + L ) 9 + w 2 

o o o o o 


]{■}-• 


substitutions 

(14*) 



(15') 


(14") 


(15") 
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Equations (14") and (15") can be solved simultaneously for 0 O and Y 0 . 


0o = 


Ki 


S 

o 





(16) 


Y 


o 



F 



(17) 


where 


K l = !o + - — 


Rewriting equation (10) 


ft 



0 




\ 

r \ 


Xi 

i 


I 

+ 

o 

1 

I ’> 


i 

4 

> 


+ bf 



(18) 


Substituting for Y 0 and 0 O in equation (18) from equations (16) 
and (17) yields a set of equations in matrix form which can be used to 
determine the total mode shapes, slopes of mode shapes, slopes of bend- 
ing mode shapes, and the natural bending frequencies. 



( 19 ) 
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The product of the row times the column matrices yields two 


square matrices which can be added to 


resulting in 


where 


{•}■*[’]{•}■ 

D 1 is the dynamic matrix 


( 20 ) 


for a free-free beam experiencing 


bending vibrations. 


The procedure to be used in the development of the influence 
coefficient matrices for digital computation can be found in Appendix A. 

An iteration procedure for obtaining higher modes is given in 
Appendix C. 


Free-Free Torsion Equations 

A set of equations for determining torsional frequencies and mode 
shapes may be written similar to those for free- free bending. The 
general equation for the angle of twist T^, of any section i is 


T 


i 




J i T i 


where 


(i - 1, 2, ... n) 



is the natural torsional frequency 

is the torsional influence coefficients, and 

is the polar mass moment of inertia. 


( 21 ) 
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Assuming that the cantilever coefficients can be determined, n 
equations may be written in matrix form as follows: 


Tic } » ^ ] [V] { T i( 


Next, the clamped end of the beam is released similar to the 
free-free bending case so that the angle of twist, free-free, may be 
written as the angle of twist, cantilever, plus some angle of twist, 
T 0 , resulting from the releasing of the clamped end. 


T 1F ■ T lc + V 


Solving for T ic above and substituting into equation (22) gives 


L iF f “ 


T 0 + “f [ *Ti ] [X s] { T lc } • 


For free-free vibrations, the following equation holds: 


J. T = 0. 
l IF 


Or, in matrix notation. 


J i 1 t if h °- 


Multiplying equation (24) by yields 
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Solving for T Q and substituting into equation (24) gives 


(28) 


where 


' I J i- 


Equation (28) in simplified form is 


{ T iF } *4 


_KI’ 


(29) 


where 




is the identity matrix, and 




" r Tc 


IS 1 

d T 


_ _ 

1 ! 

l -si 


In final form 


{ 






Equation (30) may be iterated on for the mode shapes and natural 
frequencies. 


( 30 ) 
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III. CONCLUSIONS AND RECOMMENDATIONS 

Table I illustrates that the influence coefficient method can be 
used with a reasonably small number of mass points compared to the 
Stodola method to obtain accurate frequencies and mode shapes. Slopes 
of mode shapes can also be obtained with this program, but were not 
included in this report since the mode shapes give an indication of the 
accuracy one could expect for the slopes. The accuracy of the mode 
shapes and frequencies is increased with an increase in the number of 
mass points. The first three mode shapes do not change appreciably, 
but the fourth mode is sensitive to mass point changes from 10 to 15. 

Table II compares the cantilever, free-free bending, and free- 
free torsional frequencies for the first four modes of a typical 
space vehicle whose mass and stiffness characteristics are shown in 
Figures 13 and 14. The El used in the influence coefficient method 
for each station was determined by averaging three values taken at 
1/4, 1/2, and 3/4 of the length of the mass segment. For El distribu- 
tions that vary radically over a particular mass segment, it is recom- 
mended that the reciprocal of an average value of the 1/EI distribution 
be used for the effective El. 

Figures 1 through 12 compare the first four normal modes for the 
free-free bending, cantilever bending, and the free-free torsion case 
as obtained by the two methods. Good agreement was obtained in the 
first three modes of the free-free bending case. The deviation in the 
fourth mode can possibly be explained by the fact that two extra 
"sweeping" processes were initiated before this mode was obtained. 

This arose from the fact that the El of a section of the nose was 
small in comparison to the section beginning at station X = -25 
(Figure 13). These intermediate modes or "tower modes," as they might 
be called for a vehicle with an extremely flexible tower on the nose 
of the vehicle, were not included in this report since the Stodola 
method did not indicate their existence. Excellent agreement was 
obtained for the cantilever mode shapes (Figures 5 through 8). Tor- 
sional modes were in good agreement through the second mode. A slight 
deviation occurred in the third mode and a considerable deviation in 
the fourth. A different method for obtaining cantilever bending and 
torsional influence coefficients could possibly increase the accuracy 
of the program. The more difficult variable to evaluate properly for 
each station appears to be the stiffness (El); therefore, it is recom- 
mended that various ways be tried to determine the proper El. 



TABLE OF RESULTS I 
(Uniform Beam) 
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Modes 


Free-Free Natural 
Bending Frequencies (L/ sec) 

1 

2 

3 

4 

- (145) MPTS. 

48.64 

120.52 

209.24 

305.11 

w IC - (10) MPTS. 
M’ 

48.40 

118.15 

200.42 

283.72 

- (15) MPTS. 

48.51 

119.20 

204.12 

292.49 

w IC - (20) MPTS. 

48.56 

119.58 

205.52 

295.97 


Modes 


% Variation 

1 

2 

3 

4 

ST. I. C. (10) 

0.49 

1.97 

4.22 

7.01 

ST. I. C. (15) 

0.27 

1.10 

2.45 

4.14 

ST. I. C. (20) 

0.16 

0.78 

1.80 

3.00 


ST. Stodola method 

I,C. — — Influence Coefficient method 
MPTS . .... Mass points. 


Maximum 

Normalized 

Deflections 

ST. 

I. C. (10) 

I. C. (15) 

I. C. (20) 

y il and Y 1R 

1.000 

1.000 

1.000 

0.999 

1.000^ 

0.999 

1.000 

0.999 

Y 2L and Y 2R 

1.000 

-1.000 

1.000 

-0.997 

1.000 

-0.997 

1.000 

-0.998 

Y 3L and Y 3R 

1.000 

1.000 

1.000 

0.985 

1.000 

0.986 

1.000 

+0.989 

Y 4l “ d Y 4 r 

1.000 

-1.000 

1.000 

-0.871 

1.000 

-0.951 

1.000 

-0.954 


Subscripts L and R denote left and right extremes of beam. 
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TABLE OF RESULTS II 
(Non-Uniform Beam) 


Stodola Method Modes 


Natural Frequencies Rad/sec 

i 

2 

3 

4 

(cantilever) 

2.37 

8.25 

18.41 

32.93 

(free- free) 

7.78 

18.85 

34.58 

51.66 

(torsion) 

36.41 

59.15 

86.20 

117.65 


Influence Coefficient Method (20) Pts Modes 


Natural Frequencies Rad/sec 

1 

2 

3 

4 

(cantilever) 

2.39 

8.40 

18.48 

32.48 

(free -free) 

7.95 

18.78 

33.93 

50.90 

(torsion) 

35.63 

58.97 

86.19 

115.04 


Percent Variation Modes 


U IC 

7. Var. = (1 - ~ — ) 100 

U ST 

1 

2 

3 

4 

(cantilever) 

0.84 

1.82 

0.38 

-1.36 

(free-free) 

2.19 

-0.42 

-1.88 

-1.47 

(torsion) 

-2.14 

-0.30 

0.00 

-2.22 
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FIG. 12 FOURTH FREE- FREE TORSION 
MODE FOR A NON-UNIFORM 
BEAM 
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APPENDIX A 

Cantilever Influence Coefficients for Bending 

Consider a beam divided into n equal segments of length 2 £ with 
the mass concentrated at the geometric center. Influence coefficients 
for a unit force may be written as follows: 



Using the moment area method, the deflection at station (1) due to 
a force at (1) is 


FB _ l 3 
Cl1 - • 

The centroid of the area under the M/EI diagram with respect to station 

(1) is 

kFB _ 21 
X l - — • 

It then follows that 
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for i > 2. 


The following general expressions now can be written: 


A (X - X ) 2 

^i- u L -br~ 


i=l 


i *i 


for j > 2 


and 


C^B 

r^- 




i=l 


E i *i 


for j > 2 


(P 

= -U- (f? + x ± - Xj) for i > j > 2 
j 




MB 

The influence coefficients Cjj denote the deflection of station i due 
to a unit moment at j, which produces a bending deflection and may be 
written as follows: 


(P = C ^ 

Cxi V 6E i 3E o v * 


M 

The distance from the centroid of the — diagram to station (1) is 

£j J. 


if. 




l + i_ 


2E X Ii 2E 0 Io 
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For i > 1 then 


<*f* + x i - x i>> 


and for j > 1 




f rfu_v i + _^_ 

LL E t I, J 2Ej Ij 




pMB 

JJl 


‘ * Vj] 


for j > 1 


<JO 

■ -dt <*7 + x i - V £ “ 1 > J 

i 

\ 

W for i < j. 


FB 


The influence coefficients 6j* are symmetrical with respect to 
( c ^* ref. 1); therefore, the following expressions may be written: 






and 
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The rotation of station i due to a unit moment at j, which produces 
bending, yields another set of influence coefficients 0®. Again using 
the moment area method of reference 3, 


fl MB , j l 
011 E x Ii + 2 


1 1 


E o ^o E i I i 


MB _ „MB 
’il 


0..T = for i > 1 


MB 

e jj - 2i 


t (e^z) ■ s fr 3 


for j > 1 


and 


MB MB . , . 4 

9 ij “ G j j f ° r 1 > j 


MB MB , . ^ , 

6 ij = 9 ji f ° r 1 < 


Next, shear deflections due to a unit force will be considered. 
These coefficients can be determined by multiplying two matrices con- 
taining the shear area at each station, Agi, arranged in the following 
manner: 



0 

0 

i 

1 

1 


^ A S1 

^S2 

i 

1 

1 





0 


^ A S2 


0 


0 


0 0 .. 



N x 2N-1 
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The influence coefficients 0£j may be obtained by setting up a 
diagonal matrix of the 1/A g £ values and multiplying by a triangular 
matrix as follows: 


1 

[\ t 1 

1 1 1 1 

1 

|\ l 1 


2G 

A Si 

\l 

0 1 1 1 

0 0 1 1 

• • 

• • 

• ' • 

2G 

a s 2 

\l 

N 


To determine the influence coefficients a diagonal matrix of 1/A g £ 

is multiplied times the transpose of the above triangular matrix. 



It is assumed that G is constant and therefore can be factored 

out of each of the matrices used to obtain C 1 ??, 9?? and C? 5 ?. 

ij iJ 
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APPENDIX B 

Torsional Influence Coefficients 

The torsional influence coefficients given in this appendix are 
written for a cantilever beam. R^j denotes the rotation of station i 
due to a unit torque applied at j. 


O 


i & U 21 

i un ©* 1/ 

u 

®j-i 6 

U 2 & 

£>J ©3 +1 l/\ 

l 

®- 

/ — x x -^ 

— x 2 — ► 

▲ 

w r ’ 1 

1 

• 


r T , *1 

11 G 0 Ip Q 


and 


»* . r t + - S i 

R 22 R 11 + G Ip i 


Therefore, 


T T T 

or Rt, = R 


+ X J I X J-A 

jj -J-l. J-l Oj.! I pJ .i 


j > 1 


R T . = R x . for i < j 
ij li 


and 


R ij = R jj f ° r i>J ’ 





32 


APPENDIX C 

Obtaining Higher Modes 

After the first mode shape and bending frequency have been 
determined, the higher modes and frequencies can be obtained as 
follows: 

1. Iterate on the dynamic matrix from the front and obtain 
a characteristic row. 


1 » ••• 


Ki, K 2 , K 3 , K,. 


is the original dynamic matrix. 


2. Normalize to the r^ unknown 


Kj_ Kj 

K * K ' K* 
r r t 


3. Form a square matrix with zeros for all elements in every 

. Ki 

row except the r ttl . Insert the rr- normalized row here. This will be 

Kr 

called the E x matrix. 

4. To obtain the new dynamic matrix D 2 for obtaining the second 
mode, perform the following operations: 

d 2 ] - [ D X ] 



where 


I 


is an identity matrix. 


This same procedure is used to obtain the next mode, etc. 
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APPENDIX D 

Shear Deflections Due to a Pure Moment 


The equation for displacements due to shear deformation is 




i|r S d 
T s x 

A G * 


where s is the shear distribution due to a unit load and i|r is a cor- 
rection factor used to obtain the proper shear area (A s = I A). For 

• ♦ 

a unit shear load, S * 1, the above equation can be written as follows: 



p k S d x 

Jo A ® G 


Next, consider a beam with a unit moment applied at some point j: 



dM 

Since — 
dx 


s and M’ dx 


s dx. 




M' dx 
A s G 


But M' = oo at X = Xj. However, the above integral can be evaluated by 
first considering M* as a unit finite impulse function as shown in the 
following figure. 
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